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The  radio  refractive  index  is  calculated  by  summing  the  refraction  due  to  the  infrared  resonances  of 
the  water  vapor  monomer.  Although  this  calculated  refraction  agrees  well  with  measured  values,  its 
temperature  dependence  departs  significantly  from  the  measured  variation  with  temperature.  The 
rotational  resonances  of  water  vapor  give  a  contribution  to  the  radio  refractive  index  that  does  not 
vary  precisely  as  e/T1 ;  theoretical  reasons  are  given  for  this  departure  from  the  classical  refraction  by 
permanent  dipoles. 


INTRODUCTION 

The  radio  refractive  index  equation  of  moist  air  is 
commonly  given  [e.g.,  Thayer,  1974]  by 

(n  -  1)  x  10#  =  Kt(PJT)Z,1  +  KJe/m:' 

+  K3(e/T2)Z~ 1  (1) 

where  n  is  refractive  index,  P,  is  partial  pressure  of 
dry  air,  e  is  partial  pressure  of  water  vapor,  T  is 
absolute  temperature,  and  Z,  and  Zw  are  small  cor¬ 
rection  factors  for  nonideal  gas  behavior.  The  ex¬ 
pressions  ( PJT)Z~ 1  and  ( e/T)Z~ 1  are  proportional 
to  the  number  densities  of  dry-air  molecules  and 
water  vapor  molecules,  respectively.  The  compress¬ 
ibility  factors  Za  and  Z„  are  given  by  Owens  [1967]. 
The  coefficients  K2,  and  K3  are  obtained  by 
fitting  data  to  this  equation.  This  traditional  fit  to 
'  terms  of  the  form  e/T  and  e/T2  corresponds  to  the 
classical  refraction  by  induced  dipoles  and  per¬ 
manent  dipoles,  respectively.  We  shall  be  concerned 
with  the  refraction  caused  by  water  vapor  and  hence 
with  the  corresponding  coefficients  K2  and  K3. 

Each  absorption  resonance  provides  a  nondisper- 
sive  refraction  at  frequencies  mucn  lower  than  the 
resonance  frequency.  The  electronic  transitions  at 
ultraviolet  frequencies  are  the  induced  dipole  type 
and  hence  provide  refraction  at  optical  frequencies 
•  that  varies  as  PJT  and  e/T.  Although  pure  vibra¬ 
tional  transitions  produce  a  radio  frequency  refrac- 
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tion  that  varies  as  ( e/T)  [Kan  Vleck,  1927],  the  infra¬ 
red  vibrational-rotational  transitions  do  not  entirely 
fit  either  classical  analogy,  and  Zhevakin  and  Naumov 
[1967a]  show  that  these  resonances  contribute  little 
to  the  refractive  index  at  radio  frequencies.  The  far- 
infrared  rotational  transitions  of  the  water  molecule 
are  the  permanent-dipole  type  and  can  be  expected 
to  provide  a  radio  refraction  that  principally  varies 
as  e/T2.  However,  Van  Vleck  [1932]  shows  that  for 
the  e/T 2  dependence  to  be  valid  it  is  vital  that  no 
permanent-dipole  rotational  transitions  involve 
energy  changes  of  the  order  of  the  thermal  energy 
Kt  T,  where  Kt  is  the  Boltzmann  constant  and  T  is 
the  absolute  temperature.  Equivalently,  there  must 
be  no  transitions  having  resonance  frequency  v,  (in 
units  cm"1)  of  the  order  of  KtT/hc,  where  h  is 
Planck’s  constant  and  c  is  the  speed  of  light  in 
vacuum.  Note  that  C  =  Kg/ he  =  0.695  cm-1  K'1. 
One  cm-1  is  30  GHz.  For  typical  atmospheric  tem¬ 
peratures  between  250  K  and  300  K,  the  frequency 
CT  lies  between  175  and  210  cm'1  (5.25  THz  and 
6.3  THz).  The  rotational  transitions  of  the  water 
molecule  certainly  include  strong  resonances  in  and 
near  this  frequency  interval.  Thus  there  is  theoretical 
expectation  that  the  refractive  index  of  water  vapor 
may  not  be  simply  the  sum  of  two  terms  behaving  as 
e/T  and  e/T1. 

Thayer  [1974]  hypothesized  that  the  radio  refrac¬ 
tive  index  equation  would  be  more  accurate  if  Kt 
were  taken  to  have  the  same  value  as  it  does  at  opti¬ 
cal  frequencies.  It  is  now  dear  that  this  is  tanta¬ 
mount  to  assuming  that  the  sum  of  refraction  contri¬ 
butions  from  all  the  infrared  water  vapor  resonances 
produces  K,  but  does  not  change  K2.  We  test  this 
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Fig.  1.  The  ratio  of  the  radio  refractivity  contributed  by  all 
resonances  having  resonance  frequency  less  than  a  value  on  the 
ordinate  to  the  radio  refractivity  contributed  by  all  infrared  water 
vapor  resonances.  Solid  curve  corresponds  to  T  =  273  K  while 
the  dashed  curve  is  for  T  -  333  K. 

hypothesis  by  calculating  the  sum  of  the  contri¬ 
butions  of  infrared  water  vapor  resonances  to  the 
refraction  at  radio  frequencies. 

CALCULATED  REFRACTIVE  INDEX 

In  principle,  the  refractive  index  can  be  calculated 
by  summing  the  refraction  contributions  of  all  ab¬ 
sorption  resonances  [Zhevakin  and  Naumov,  1967 a, 
b;  Birnbaum ,  1953].  Birnbaum  [1953]  shows  that  for 
any  sharp  absorption  resonance  the  Kramers-Kronig 
relationship  implies  a  refraction  in  the  far  wings 
having  a  frequency  dependence  given  by  n,  oc  1/ 
(v?  —  v2),  where  n,  is  the  refractive  index  caused  by 
the  resonance  labeled  by  the  index  i,  v  is  frequency, 
and  v,  is  the  resonance  frequency.  This  result  applies 
for  either  v  «  v,  or  v  »  v, .  Since  this  result  applies  to 
sharp  absorption  resonances  in  general,  it  certainly 
applies  for  any  particular  line  shape  model  having  a 
sharp  absorption  line  shape  (e.g.,  the  Gross  or  Van 
Vleck-Weisskopf  line  shapes).  Eac'.i  resonance  con¬ 
tributes  a  refraction  n,  that  is  nondispersive  in  the 
limit  of  low  electromagnetic  frequency  v-*0;  each 
line  contributes  an  amount  [Hill  and  Clifford,  1981] 

n,  -  S,Q/2n2vf  (2) 

where  S,  is  the  integrated  absorption  strength  of  the 
resonance,  v,  is  its  line  frequency,  and  Q  is  the 
number  density  of  absorbing  molecules.  Any  particu¬ 
lar  line  shape  model  that  is  sufficiently  sharp  in  the 
sense  of  Birnbaum  [1953]  (e.g.,  the  Gross  or  Van 


Vleck-Weisskopf  line  shapes)  must  give  equation  (2). 
The  strongest  (largest  S()  and  lowest  frequency 
(smallest  v,)  resonances  contribute  most.  Therefore, 
although  the  strong  water  vapor  resonances  in  the 
vicinity  of  the  frequency  CT  may  not  give  a  temper¬ 
ature  dependence  behaving  like  e/T2,  the  equally 
strong  but  lower  frequency  resonances  may  dominate 
the  radio  wave  refraction  and  approximately  restore 
this  temperature  dependence. 

The  line  strength  has  temperature  dependence  aris¬ 
ing  from  the  partition  function  and  the  difference  of 
two  Boltzmann  distributions.  For  the  water  molecule 
this  temperature  dependence  is 

S,  x  r-wfexp  (-Et/CT)  -  exp  (-fif'/CT)]  (3) 

where  E[  and  E\  are  the  values  (in  units  cm  * l)  of  the 
lower  and  upper  state  energies  responsible  for  reson¬ 
ance  i  and  C  =  0.695  cm " 1  K~\ 

We  have  calculated  the  radio  frequency  refraction 
of  the  infrared  resonances  of  water  vapor  by  sum¬ 
ming  the  contributions  of  all  water  vapor  lines  in  the 
Air  Force  Geophysics  Laboratory  (AFGL)  line  pa¬ 
rameter  compilation.  This  compilation  extends  from 
lines  at  radio  frequencies  to  those  in  the  visible.  All 
water  vapor  lines  are  included  in  the  sum;  no  lines 
have  been  excluded,  no  matter  how  weak  or  how 
high  in  line  frequency.  The  calculation  was  per¬ 
formed  for  several  humidities  from  0%  to  100%  for 
each  of  many  temperatures  from  -  10°C  to  63°C  and 
1 10°C  to  250°C  and  for  total  pressures  of  0.2  atm,  0.8 
atm,  and  1.0  atm.  As  expected,  for  fixed  temperature 
and  humidity  the  refraction  by  water  vapor  was 
found  to  be  independent  of  total  pressure.  Also  as 
expected,  the  calculated  refraction  was  proportional 
to  the  humidity.  The  Gross  line  shape  having  humi¬ 
dity  and  pressure-dependent  line  width  was  used,  as 
by  Hill  et  al.  [1980].  This  line  shape  gives  (2)  in  the 
low-frequency  limit.  A  1-m  radio  wavelength  was  as- ' 
sumed;  for  such  a  low  frequency  there  is  no  practical 
difference  between  the  complete  Gross  formula  and 
(2). 

Suppose  we  begin  summing  the  refraction  with  the 
lowest  frequency  in  the  line  parameter  compilation 
and  proceed  to  include  higher  and  higher  frequency 
resonances.  How  does  the  refraction  accrue  as  more 
resonances  are  included,  and  how  high  in  frequency 
must  we  sum  before  the  total  refraction  is  ap¬ 
proached?  Figure  1  illustrates  the  answer  by  showing 
the  ratio  of  the  partial  sum  to  the  total  sum.  Al¬ 
though  the  graph  only  shows  frequencies  less  than 
400  cm'1  (i.e.,  12  THz),  the  total  sum  includes  all 
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lines  from  radio  to  visible  frequencies.  In  principle, 
the  curves  in  Figure  1  should  step  upward  at  the 
resonance  frequency  of  each  line;  in  fact,  the  curves 
are  smoothed  by  plotting  point  to  point  at  intervals 
of  10  cm" 1 .  Consider  the  proof  by  Van  Vleck  [1932] 
that  rotational  resonances  should  sum  to  provide  a 
radio  refraction  that  varies  as  (e/7’2);  this  proof  re¬ 
quires  neglecting  Q.S^JCT)1  relative  to  unity.  If 
0 .5(vJCT)2  is  less  than  0.1,  then  v,  <  85  cm-1  for 
T  =  273  K.  Figure  1  shows  that  less  than  half  of  the 
radio  refraction  arises  from  rotational  resonances  sa¬ 
tisfying  O.S^vJCT)2  <  0.1.  Thus  we  should  not  be  sur- 
r  ts--i  if  the  radio  refraction  caused  by  rotational 
resonances  does  not  vary  precisely  as  (e/7’2). 

In  a  sense,  a  measurement  of  radio  refractive  index 
is  an  empirical  sum  of  the  refraction  contributions  of 
all  the  infrared  and  ultraviolet  resonances.  It  is  natu¬ 
ral  to  treat  our  calculated  refractive  indexes  as  a  si¬ 
mulated  set  of  experimental  data.  Thus  we  do  as  the 
experimentalists  have  done;  that  is,  we  fit  the  calcu¬ 
lated  refraction  to  the  form  K2(e/T)Z~ 1  + 
K3{e/T2)Zw  1  as  suggested  by  (1).  This  functional 
form  fits  the  refraction  very  well;  the  rms  error  was 
only  a  few  thousandths  of  one  percent.  The  temper¬ 
ature  range  0°-60°C  corresponds  to  that  used  in  the 
experiments  by  Boudouris  [1963];  fitting  to  this  tem¬ 
perature  range  gave  the  following  infrared  contri¬ 
butions  to  K2  and  K3  : 

K2  -  33.30  ±  1  K/mbar  (4) 

K,  =  (3.5833  ±  0.0032)  x  105  KJ/mbar  (5) 

Equations  (4)  and  (5)  result  from  the  sum  of  all  infra¬ 
red  resonances  in  the  AFGL  compilation.  Note  that 
the  total  value  of  K2  at  radio  frequencies  is  the  sum 
of  the  above  infrared  contribution,  33.30  K/mbar, 
and  the  extrapolation  to  infinite  wavelength  of  the 
empirically  known  value  of  K2  at  optical  frequencies, 
64.79  ±0.08  K/mbar  IThayer,  1974];  the  latter  is 
caused  by  the  ultraviolet  electronic  resonances  of  the 
water  molecule. 

We  examined  fitting  the  calculated  refraction  to 
several  other  temperature  ranges.  For  the  temper¬ 
ature  range  1 10°-250°C,  which  roughly  corresponds 
to  the  higher  temperature  range  of  experiments  re¬ 
viewed  by  Bimbaum  and  Chatterjee  [1952],  we  found 
that  the  infrared  contribution  to  K2  is  2%  larger 
than  (4)  and  K3  is  0.1%  smaller  than  (5).  For  the 
temperature  range  - 10°  to  40°C,  the  infrared  contri¬ 
bution  to  AC j  is  2%  larger  than  (4),  and  AC3  is  0.05% 
smaller  than  (5).  These  deviations  are  within  the  error 
given  in  (4)  and  (5). 
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The  choice  of  the  estimated  errors  in  (4)  and  (5) 
deserves  explanation.  Unlike  experiments  which  have 
random  variation,  we  are  fitting  to  a  smooth  function 
of  temperature.  We  must  choose  a  nonstatistical 
measure  of  the  error  in  K2  and  AC3 .  Our  choice  is  as 
follows:  We  first  find  the  best  fit  that  the  computer 
can  generate;  then  we  make  a  small  change  in  K2 
and  refit  by  changing  only  Af3,  but  not  K2.  We  con¬ 
tinue  to  change  K2  and  refit  K3  until  the  rms  error  in 
the  refraction  is  twice  the  rms  error  of  the  best  fit. 
The  errors  stated  in  (4)  and  (5)  are  the  difference 
between  the  values  of  K2  and  AC3  from  the  best  fit 
and  those  that  double  the  rms  error  in  the  refraction. 
This  definition  of  error  gives  a  conservative  estimate 
that  accounts  for  the  fact  that  a  decrease  in  K2  can 
be  approximately  compensated  by  an  increase  in  K3 
and  vice  versa.  If  all  the  rms  error  in  the  fit  to  the 
refraction  were  blamed  on  only  one  of  K2  or  K3 ,  but 
not  the  other,  then  the  errors  would  be  only  about 
one-tenth  the  errors  stated  in  (4)  and  (5).  The  differ¬ 
ence  between  minimizing  rms  error  as  opposed  to 
minimizing  percent  rms  error  changes  K2  and  K3  by 
about  one-tenth  the  error  stated  in  (4)  and  (5).  The 
error  in  the  fitting  is  not  caused  by  computer  round¬ 
off;  it  is  caused  by  the  computed  refraction  having 
temperature  dependence  that  is  not  adequately  de¬ 
scribed  by  (1). 

The  line  strengths  are  given  to  three  decimal  places 
in  the  line  parameter  compilation.  By  allowing  in¬ 
creases  in  K3  to  partially  compensate  for  decreases  in 
K2  (and  vice  versa),  three  significant  figures  in  the 
calculated  refraction  imply  errors  in  K2  and  K3  that 
are  about  10  times  larger  than  those  stated  in  (4)  and 
(5).  As  will  be  seen  in  the  discussion  section,  such 
large  errors  equal  the  experimental  errors  in  K2  and 
K3 .  Whether  or  not  the  three  decimal  places  are  all 
accurate  is  another  question  entirely. 

We  evaluate  the  relative  contributions  to  K2  and 
K3  from  the  rotational  resonances  versus  the 
vibrational-rotational  resonances  by  separately  cal¬ 
culating  the  refraction  from  all  infrared  resonances  at 
frequencies  in  excess  of  1000 cm-1  (corresponding  to 
a  10- jim  wavelength)  and  separately  for  all  reson¬ 
ances  at  frequencies  less  than  1000  cm'1.  The  arbi¬ 
trary  boundary  of  1000  cm'1,  which  lies  in  the 
10-^m  window,  suffices  to  separate  the  radio  refrac¬ 
tion  caused  by  the  vibrational-rotational  transitions 
at  higher  frequencies  from  the  pure  rotational  transi¬ 
tions  at  lower  frequencies.  We  find  that  nearly  all  of 
both  K2  and  K3  is  caused  by  rotational  resonances. 
The  vibrational-rotational  resonances  produce  only 
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TABLE  1.  Values  of  Kj  and  K  ,  for  Radio  Frequencies 


Boudouris 

72  ±  10.5 

(3.754  ±  0.03)  x  10* 

Birnbaum  and  Chatterjee 

71.4  ±  5.8 

(3.747  ±  0.029)  x  10* 

Thayer 

64.79  ±  0.08 

(3.776  ±  0.004)  x  10* 

Calculated 

98  ±  1 

(3.583  ±  0.003)  x  10* 

Values  of  K i  are  in  degrees  Kelvin  per  millibar;  values  of  K , 
are  in  degrees  Kelvin  squared  per  millibar. 


5.8%  of  the  infrared  contribution  to  K2  and  decrease 
K3  by  only  0.0009%.  Specifically,  for  the  fit  to  the 
temperature  range  0°-60°C,  the  infrared  water  vapor 
resonances  at  v,  >  1000  cm-1  give  a  contribution  to 
K2  of  1.92  ±  0.02  K/mbar  and  a  contribution  to  K3 
equal  to  —  30.9  ±  5  K2/mbar.  Notice  that  the  net 
refraction  by  these  resonances  is  positive.  We  have 
no  objection  to  a  negative  contribution  to  K3  be¬ 
cause  we  do  not  force  a  classical  interpretation  on 
the  refraction  by  the  vibrational-rotational  reson¬ 
ances. 

We  have  investigated  forcing  the  rotational  reson¬ 
ances  to  contribute  to  only  one  term  in  (1).  For  the 
temperature  range  0°-60°C,  we  fitted  the  calculated 
refraction  from  all  resonances  having  v(  <  1000  cm  ~ 1 
to  the  form  K'3  (e/Tr)  and  obtained 

K’3  m  (3.1748  ±  0.0083)  x  10s  K'/mbar  (6) 

p  «=  1.97421  ±  0.00046  (7) 

The  errors  are  defined  analogous  to  the  errors  stated 
in  (4)  and  (5).  This  one-term  fit  gave  an  rms  error  of  a 
few  thousandths  of  one  percent  in  the  refraction;  this 
is  roughly  the  same  error  as  obtained  by  fitting  with 
the  two  terms  in  (1).  We  fitted  the  calculated  refrac¬ 
tion  from  all  resonances  v,  >  1000  cm-1  to  the  form 
K’t (e/T*)  and  obtained 

K't  m  1.319  ±  0.038  K'/mbar  (8) 

q  «  1.0565  ±  0.0050  (9) 

This  two-parameter,  one-term  fit  gave  a  0.02%  rms 
error  in  the  refraction  caused  by  those  resonances 
having  v,  >  1000  cm-1.  In  comparison,  the  two- 
parameter,  two-term  fit  in  (1)  gave  0.03%  rms  error. 

DISCUSSION 

Table  1  gives  the  values  of  K2  and  K3  resulting 
from  measurements  by  Boudouris  [1963],  the  review 
of  measurements,  including  their  own,  by  Bimbaum 
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and  Chatterjee  [1952],  the  values  reported  by  Thayer 
[1974],  who  assumed  K2  equals  its  value  at  optical 
frequencies,  and  the  calculations  described  in  the  pre¬ 
vious  section.  Birnbaum  and  Chatterjee  [1952]  col¬ 
lected  the  K2  and  K3  values  from  several  experi¬ 
ments;  they  averaged  the  K2  values  weighted  in¬ 
versely  proportional  to  the  square  of  the  probable 
errors  and  averaged  the  K3  values  by  the  same 
method.  The  values  of  K2  and  K3  given  by  Smith  and 
Weintraub  [1953]  are  based  on  those  by  Bimbaum 
and  Chatterjee  [1952].  On  the  basis  of  their  recent 
measurement  of  water  vapor  refraction,  Liebe  et  at. 
[1977]  give  values  of  K2  and  K3  that  are  nearly 
identical  to  those  of  Smith  and  Weintraub  [1953]. 
The  appendix  discusses  the  significantly  different 
values  of  K2  and  K3  obtained  by  Hasegawa  and 
Stokesberry  [1975].  The  experimental  values  of  K2 
and  K3  seem  consistent  and  well-represented  by  the 
upper  two  entries  in  Table  1. 

Now  compare  the  experimental  values  in  Table  1 
(those  of  Boudouris  and  of  Birnbaum  and  Chattetjee) 
with  the  calculated  values.  Similar  to  the  calculations 
by  Zhevakin  and  Naumov  [1967a],  our  calculated  re¬ 
fractive  index  is  within  a  few  percent  of  measured 
values.  At  -  20°C,  the  refractivity  from  our  calcu¬ 
lated  values  of  K2  and  K3  is  2.65%  less  than  that 
obtained  using  the  K2  and  K3  values  of  Boudouris. 
This  underestimate  decreases  with  increasing  temper¬ 
ature  to  0.76%  at  260°C.  At  20°C  our  calculated  K2 
and  Kj  values  give  2.4%  less  refractivity  than  K2 
and  K3  from  Boudouris.  For  20°C,  Table  2  by 
Boudouris  [1963]  summarizes  experimental  error  in 
refractivity  for  nine  experiments;  these  errors  average 
to  1.2%.  This  is  one-half  the  discrepancy  between 
calculated  and  measured  refractivity  at  20°C.  How¬ 
ever,  the  apportionment  between  the  coefficients  K2 
and  K3  differs  from  the  measured  results  by  several 
standard  deviations.  For  instance,  assuming  that  K2 
has  an  average  experimental  standard  deviation  of  8 
(see  Table  1),  then  the  experimental  value  of  K2  is  3.2 
standard  deviations  smaller  than  its  calculated  value. 
Likewise,  the  experimental  value  of  K3  is  5.7  stan¬ 
dard  deviations  larger  than  its  calculated  value.  Thus 
if  one  plots  (n  -  i)/(e/T)  as  a  function  of  T"  \  then 
the  calculated  curve  (a  straight  line  to  excellent  ap¬ 
proximation)  is  significantly  different  from  the  mea¬ 
sured  curve  because  the  slope  K3  and  intercept  Kx  of 
the  straight  line  fits  differ  by  many  standard  devi¬ 
ations.  This  is  a  genuine  disagreement  betweeu 
theory  and  experiment 

One  possible  explanation  for  the  discrepancy  is 
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that  the  line  parameter  compilation  is  incorrect,  in 
that  Si/vf  is  inaccurate  for  some  lines.  Considering 
our  previous  comments  (see  the  introduction)  on  Van 
V leek's  [1932]  derivation,  the  rotational  resonances 
most  likely  to  provide  a  deviation  from  the  K3(e/T2) 
dependence  are  those  having  frequencies  near  v,  = 
KgT/hcz:  190  cm-1  (i.e.,  5.7  THz).  Since  the  term 
K2(e/T)  is  much  smaller  than  the  K3(e/T2)  term, 
even  a  small  deviation  from  a  1/T 2  dependence  by 
the  rotational  contribution  will  cause  a  large  contri¬ 
bution  to  K2. 

A  second  possible  explanation  of  the  discrepancy  is 
that  water  clusters  contribute  to  the  measured  refrac¬ 
tion,  whereas  we  have  only  calculated  the  monomer 
contribution.  In  principle,  it  is  possible  to  calculate 
the  refraction  contributed  by  a  model  absorption 
spectrum  of  the  water  dimer  (the  dimer  absorption 
spectrum  is  given  by  Bohlander  [1979]).  However,  it 
is  not  possible  to  calculate  the  refraction  contributed 
by  the  broadband  empirical  excess  absorption,  for 
two  reasons.  First,  the  Kramers-Kronig  relation  is 
sensitive  to  the  way  in  which  this  absorption  de¬ 
creases  at  high  frequencies.  Second,  the  excess  ab¬ 
sorption  is  defined  as  the  measured  absorption  less 
the  monomer  absorption  predicted  with  a  particular 
line  shape.  To  the  extent  that  the  tails  of  this  particu¬ 
lar  monomer  line  shape  are  inaccurate,  the  excess 
absorption  cannot  be  entirely  attributed  to  water 
clusters;  but  one  must  accurately  remove  monomer 
absorption  from  the  excess  absorption  because  the 
refraction  contribution  of  the  monomer  is  fully  ac¬ 
counted  for  in  our  calculations. 

The  discrepancy  is  not  caused  by  our  calculation 
using  a  1-m  radio  wavelength  (300  MHz)  as  opposed 
to  the  higher  radio  frequencies  typically  used  in  ex¬ 
periments  (see  Table  2  by  Boudouris  [1963]).  The 
dispersion  in  (n  -  1)  x  106  and  Ar  is  shown  graphi¬ 
cally  by  Hill  and  Clifford  [1981],  where  AT  s 
-K3(e/T2);  thus  the  dispersion  in  AT  is  a  measure  of 
that  in  K3 .  The  dispersion  in  K2  is  easily  calculated 
from  that  in  (n  -  1)  x  106  and  X3 .  From  zero  fre¬ 
quency  to  60  GHz  the  graphs  by  Hill  and  Clifford 
[1981]  imply  that  K3  increases  by  0.0096  x  10s 
K2/mbar  and  K2  decreases  by  1.6  K/mbar.  These 
changes  in  K2  and  K3  are  considerably  smaller  than 
their  experimental  uncertainties. 

Now  consider  Thayer's  values  of  K2  and  K3  in 
Table  1.  Those  values  are  based  on  his  hypothesis 
that  K2  should  be  the  same  at  radio  wavelengths  as 
at  visible  wavelengths,  namely  K2  =  64.79  K/mbar. 
We  have  seen,  however,  that  the  infrared  resonances 


of  water  vapor  are  theoretically  expected  to  give  a 
contribution  to  K2 ;  and  the  calculated  values  of  K2 
show  that  the  infrared  resonances  of  water  vapor  can 
give  a  significant  contribution  to  K2  (namely  a  calcu¬ 
lated  contribution  of  33.3  K/mbar  in  (4)).  Thus 
Thayer’s  hypothesis  is  unfounded.  We  suggest  that 
one  should  use  the  entirely  empirical  coefficients  (e.g^ 
those  of  Boudouris  [1963]).  Of  course,  the  term 
K2{e/T)  is  only  needed  for  precise  work.  For  most 
practical  applications,  this  term  may  have  its  ap¬ 
proximate  value  introduced  into  the  K3{e/T2)  term 
to  yield  the  simpler  radio  refractive  index  equations 
discussed  by  Smith  and  Weintraub  [1953]  or  Bean 
and  Dutton  [1966], 

We  conclude,  quite  opposite  to  the  assumption  by 
Zhevakin  and  Naumov  [1967a],  that  the  infrared  con¬ 
tribution  to  K  2  is  dominated  by  the  rotational  reso¬ 
nances,  not  by  the  vibrational-rotational  resonances. 
We  noted  in  the  introduction  that  the  rotational  res¬ 
onances  of  water  vapor  need  not  provide  a  temper¬ 
ature  dependence  proportional  to  1/T2.  In  fact,  Zhe¬ 
vakin  and  Naumov  [1967a]  calculated  the  temper¬ 
ature  dependence  of  radio  refractive  index  caused  by 
the  rotational  resonances;  they  obtained  1/T1  97.  If 
one  fits  1/T197  to  the  form  (K2/T)  +  (K3/T2),  one 
obtains  a  value  of  K2  that  is  roughly  the  magnitude 
of  the  K2  contribution  that  we  found  from  the 
rotational  resonances.  Thus  the  contribution  of  the 
rotational  resonances  to  K2  is  very  sensitive  to  the 
deviation  from  1/T2  of  their  refractive  temperature 
dependence. 

One  could  force  the  rotational  transitions  to  con¬ 
tribute  to  only  one  term  in  the  radio  refractive  index 
equation  by  fitting  their  radio  refraction  to  K’3(e/T'). 
Then  p  would  be  close  to,  but  not  equal  to,  2  (just  as 
Zhevakin  and  Naumov  [1967a]  found  p  ~  1.97).  Our 
fit  to  this  functional  form  gives  the  K'3  and  p  values 
in  (6)  and  (7).  One  would  then  likewise  require  that 
the  vibrational-rotational  transitions  contribute  to 
only  one  term  by  fitting  their  radio  refraction  contri¬ 
bution  to  K'2{e/T*).  Then  q  would  be  nearly,  but  not 
precisely,  unity.  Our  fit  gives  the  K'2  and  q  values  in 
(8)  and  (9).  We  must  view  this  type  of  fit  to  be  in  no 
fundamental  sense  superior  to  that  in  (IX  and  this  fit 
obscures  comparison  of  the  calculated  and  measured 
refractive  indexes.  Since  the  rotational  transitions  do 
not  provide  K3  exclusively  and  the  vibrational- 
rotational  transitions  do  not  contribute  to  K2  exclu¬ 
sively,  we  must  view  (1)  as  a  tow-order  polynomial  fit 
in  inverse  powers  of  T  to  a  general  function  of  tem¬ 
perature. 
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In  fact,  Van  Vleck  [1927]  shows  that  the  quantum 
mechanical  correction  to  (1)  has  the  form  K^fe/T*) 
+  Kf{e/T 4)  +  •  •  ■  .By  including  the  term  K4(e/T*) 
in  (1)  and  using  Van  Vleck’s  expression  for  KA,Birn- 
baum  and  Chatterjee  [1952]  (in  their  appendix)  show 
that  the  value  of  K2  is  reduced  from  its  value  de¬ 
duced  using  (1)  without  the  correction  term.  They 
show  that  this  reduced  value  of  K2  is  less  than  the 
extrapolation  to  zero  frequency  of  the  optical- 
frequency  value  of  K2;  this  is  impossible  since  the 
(positive)  infrared  contribution  to  K2  must  be  added 
to  its  extrapolated  optical  value  to  obtain  the  radio 
frequency  value  of  K2.  Since  Van  Vleck's  derivation 
assumes  that  there  are  no  strong  resonances  having 
energy  changes  of  the  order  of  the  thermal  energy 
KBT  and  since  we  have  seen  in  the  introduction  that 
the  water  molecule  violates  this  assumption,  it  fol¬ 
lows  that  Van  Vleck’s  value  for  Kt  should  not  be 
applied  to  water  vapor  and  the  discrepancy  noted  by 
Birnbaum  and  Chatterjee  [1952]  is  removed. 

Summing  the  refraction  by  absorption  resonances 
can  shed  light  on  the  radio  refraction  of  dry  air  as 
well  as  water  vapor.  Liebe  [1981]  calculates  refrac¬ 
tion  by  oxygen  resonances  having  resonance  fre¬ 
quency  less  that  1000  GHz.  These  magnetic  dipole 
transitions  of  molecular  oxygen  are  expected  to  be 
the  main  contributors  to  the  magnetic  susceptibility 
of  air  at  radio  frequencies.  Liebe’s  calculated  refrac¬ 
tion  caused  by  these  oxygen  resonances  accounts  for 
three  fourths  of  the  observed  magnetic  susceptibility 
of  air;  reasons  for  the  remaining  difference  in  calcu¬ 
lated  versus  observed  susceptibility  are  unknown. 

APPENDIX 

Hasegawa  and  Stokesberry  [1975]  obtained  values 
of  K2  and  K}  that  are  significantly  different  from 
those  of  Boudouris  [1963]  and  Birnbaum  and 
Chatterjee  [1952].  Hasegawa  and  Stokesberry  [1975] 
used  the  Bimbaum-Chatterjee  method  of  averaging 
K2  and  K2  from  several  different  experiments.  Unlike 
Birnbaum  and  Chatterjee,  Hasegawa  and  Stokes¬ 
berry  included  experiments  that  covered  only  very 
limited  ranges  in  temperature,  namely  the  two  experi¬ 
ments  by  Essen  and  Froome  [1951]  and  Essen 
[1953].  Because  of  their  limited  temperature  range, 
Essen  and  Froome  [1951]  assumed  K2  equal  to  its 
value  at  optical  frequencies  (an  inaccurate  assump¬ 
tion)  and  fitted  only  K3  to  their  data.  Because  of 
their  small  probable  error,  these  two  experiments 


overwhelmingly  dominated  the  average  K2  and  K3 
determined  by  Hasegawa  and  Stokesberry  [1975].  Be¬ 
cause  of  the  very  limited  temperature  range  of  these 
two  experiments  and  the  inappropriate  assumption 
regarding  K2,  we  feel  that  the  average  values  of  K2 
and  K3  of  Hasegawa  and  Stokesberry  are  not  accu¬ 
rate.  It  seems  preferable,  but  perhaps  impossible,  to 
fit  (1)  to  a  compilation  of  available  experimental 
values  of  refractivity,  rather  than  to  average  the  coef¬ 
ficients  of  fits  made  to  each  individual  data  set. 
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